Numerical study of the frustrated ferromagnetic spin-^ chain 
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The ground state phase diagram of the frustrated ferromagnetic spin-1/2 chain is investigated 
using the exact diagonalization technique. It is shown that there is a jump in the spontaneous 
magnetization and the ground state of the system undergos to a phase transition from a ferromag- 
netic phase to a phase with dimer ordering between next-nearest-neighbor spins. Near the quantum 
transition point, the critical behavior of the ground state energy is analyzed numerically. Using 
a practical finite-size scaling approach, the critical exponent of the ground state energy is com- 
puted. Our numerical results are in good agreement with the results obtained by other theoretical 
approaches. 

PACS numbers: 75.10.Jm, 75.10.Pq 
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I. INTRODUCTION 

The physics of frustrated quantum spin systems have 
been attracted much interest from experimental and the- 
oretical points of view. The spin-^ Hamiltonian of the 
frustrated model on a periodic chain of A^-sites is 



N 



H — '^^{Ji'§ n.~§ n+l + J2~§ n-~§ n+2)i 



(1) 



where ~§ n represents the 5* = ^ operator at the n-th 
site, and Ji, J2 are the nearest- neighbor (NN) and next- 
nearest- neighbor (NNN) interactions. We introduce the 
parameter a = for convenience. 

This model with NN and NNN antiferromagnetic in- 
teractions (Ji, J2 > 0) is well studiedii^i^iii^i^ii^. This 
chain is well known to display a quantum phase transi- 
tion from a gapless, translationally invariant state with 
algebraic spin correlations (the spin fluid phase) to dimer 
gapful state at ac — 0.2411—. At the Majumdar-Ghosh 
points, i.e. ai a — 0.5, the ground state is exactly solv- 
able. It is a doubly degenerate dimer product of singlet 
pairs on neighboring sites. In general, the ground state is 
doubly degenerate for a > Uc. For large J2 (a > 0.5) an 
incommensurate phase appears in the ground state phase 
diagram^'^. The behavior of frustrated chains in the pres- 
ence of a uniform magnetic field was first studied by R. 
ChitrEfc^. Recently the effect of a uniform magnetic field 
on the Ji — J2 model has been discussed^°. They have 
shown that a chiral phase emerges in isotropic frustrated 
spin chains as well, if they are subject to a strong exter- 
nal magnetic field. When Ji > and J2 < (AF-F) , the 
system is believed to be in a gapless antiferromagnetic 
phase for any permissible values of Ji and J2. 

Relatively little attention has been paid to frustrated 
ferromagnetic chains, i.e., Ji < and J2 > 0. From 
experimental point of view, the recent discovery of 
materials are described by parameters with this com- 
binations of signs. The Rh2Cu2Mo-iOi2 is believed 
to be describediiii^ by Ji ~ — 3J2, and LiCuVo4^ 
which lies in a different parameter regime with Ji ~ 
— 0.3J2^^. A recent studj^ of the thermodynamics 



of the model (1) was motivated by the experimen- 
tal results for Rh2Cu2Mo^Oi2. From theoretical point 
of view, the later model has been subject of many 
studies^ii^ii^iiiii^ii^iSa. The complete picture of the 
phases of this model as a function of the frustration pa- 
rameter a is unclear up to now. 

In the case of Ji < and J2 > (F-AF) with 
< a < -J, the ground state is fully ferromagnetic 
and hes in the subspace Stot — N/2 with the degener- 
acy -I- 1, and becomesi^ an [S = 0) incommensurate 
singlet state^iiS^ for a > also the lattice translational 
symmetry is thought to be broken. It is suggested that 
in this incommensurate singlet state, the gap is strongly 
suppresed^^. At the critical point ac = j, two distinct 
configurations with the energy 
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N I Ji 



(2) 



i is known 



are the ground states^'*. One is fully ferromagnetic with 
Stot — N/2, the other is a singlet state with Stot — 0. 
The wave function of the singlet state at a, 
exactly^^^ . 

In the vicinity of the critical point Ofc = ^ , at < 7 <C 
1 (7 = a — -j) the singlet ground state energy behaves as 
Eq ^7'^, where /3 is a critical exponent. The classical ap- 
proximation gives P = 2. The spin- wave theory as well as 
some other approximations®'^^ do not change this criti- 
cal exponent. In Ref. [26], using variational approaches, it 
has been shown that the quantum fluctuations definitely 
change the classical critical exponent. They conjectured 
that strong quantum fluctuations change the critical ex- 
ponent and /? = |. In a recent work, Dmitriev et.al.— 
have studied the properties of this model using the per- 
turbation theory (PT) in the small parameter charac- 
terizing the deviation from the transition point. They 
considered the Hamiltonian (1) as 

H = Ho + V-, 

Hq = — ~^ n-~^ 71+1 + 7 ~^ n-~^ n+2 



(3) 
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with a small parameter < 7 ^ 1. Since the pertur- 
bation Vj conserves the total spin S"^, the PT to the 
lowest singlet state | ip) of the Hamiltonian Hq involves 
only singlet excited states. They showed that the PT al- 
low them to estimate the critical exponent of the ground 
state energy as 



^0(7)^ -^7" /3 = 5/3, 



(4) 



which is in good agreement with their previous result^^. 
On the other hand, they have claimed that the exact 
diagonalization of finite chains shows a complicated ireg- 
ular size dependence of the ground state energy, which 
makes the numerical estimation of the critical expo- 
nent f3 impossible?^. In a very recent work, the ground 
state phase diagram of the spin-1/2 zigzag chain with 
weakly anisotropic ferromagnetic NN and antiferromag- 
netic NNN interactions is studied—. It is shown that 
the ground state phase diagram consists of the fully po- 
larized ferromagnetic, the commensurate spin-liquid and 
the incommensurate phases. 

In this paper, we present our numerical results on the 
ground state phase diagram of the ID frustrated ferro- 
magnetic spin-1/2 model. Our results are obtained us- 
ing the exact diagonalization technique. In section 2, we 
present the results of exact diagonalization calculations 
on the ground state phase diagram of the model. In sec- 
tion 3, we discuss a practical finite-size scaling approach 
and find the critical exponent of the ground state energy 
in the vicinity of the critical point ac = \- Finally, the 
summary and conclusions are presented in section 4. 



II. THE GROUND STATE PHASE DIAGRAM 
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FIG. 1: The spontaneous magnetization as a function of 
the parameter a for different chain lengths A'^ = 26, 28, 30. 



To recognize the different phases induced by the NNN 
exchange interaction in the ground state phase diagram, 
we have implemented the Lanczos algorithm of finite size 
chains to calculate the order parameters and the various 
spin correlation functions. The first insight into the na- 
ture of the different phases can be obtained by studying 
the uniform magnetization 
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x,y,z 



(5) 



An inportant goal in the study of quantum spin sys- 
tems is the search for novel states emerging from com- 
peting interactions in the ground state phase diagram. 
In particular, the study of continuous phase transitions, 
has been one of the most fertile branches of theoretical 
physics in the last decades. Each phases can usually be 
characterized by an order parameter. Often, the choice 
of an order parameter is obvious, however in some cases 
finding an appropriate order parameter is complicated. 
As we mensioned, the complete picture of the phases of 
this model as a function of the frustration parameter a is 
not completely clear. It is known that the ground state 
is ferromagnetic at < a < ac and a second order phase 
transition happens to the incommensurate singlet phase. 

In order to explore the nature of the spectrum and 
the phase transition, we used the Lanczos method to di- 
agonalize numerically finite (up to A'^ = 30 sites) chain 
systems. The energies of the few lowest eigenstates 
were obtained for chains with periodic boundary con- 
ditions. The Lanczos method and the related recursion 
methods,— i^i^i^ possibly with appropriate implemen- 
tations, have emerged as one of the most important com- 
putational procedures, mainly when a few extreme eigen- 
values are desired. 



where the notation (...) represent the expectation value 
at the lowest energy state. 

In Fig. [1] we have plotted the spontaneous magne- 
tization, vs a for the chain of different lengths 
N — 26, 28, 30. For arriving at this plot we consid- 
ered I Ji 1= 1 and different values of the parameter 
< a < 0.5. One of the most interesting known prop- 
erties of this model is that the magnetization as a func- 
tion of applied magnetic field displays a jump for certain 
parameteraiSi^i^^. It can be seen that the spontaneous 
magnetization , remains close to the saturation value 
for < a < ttc. This behavior is in agreement with ex- 
pectations based on the general statement that for values 
of the parameter < a < ac, the ground state is in the 
gapped ferromagnetic phase. At a critical value a = ac, 
the spontaneous magnetization jumps to zero. However, 
we observe that the metamagnetic phase transition oc- 
curs also in the absence of the external uniform magnetic 
field. The zero value of the spontaneous magnetization 
in the region a > ac, shows that the ground state of the 
model is not magnetic. 

To disply the ground state magnetic phase dagram of 
the model we have calculated the dimer order parame- 
ters. Because of two type of coupling constants, we in- 
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FIG. 2: (a) The F-dimer order parameter dp as a function of 
parameter a for different chain lenghts A'^ — 20, 26. (b) The 
AF-dimer order parameter cIaf as a function of parameter a 
for different chain lenghts TV = 20, 26. 



troduce two kind of dimerization as 

d-F = -^Yl {~^n-^,i+i) , (6) 



a > ac and enough large values of parameter a, F-dimer 
order parameter is slightly more than zero dp ~ 0.05 
but the AF-dimer order parameter is less than saturation 
value (—0.75) cIaf ^ —0.45. Thus, increasing the antifer- 
romagnetic exchang J2 from critical value ac, quantum 
fluctuations suppress the ferromagnetic ordering and the 
system smoothly undergoes a transition from a ferromag- 
netic phase into a phase with the dimer ordering between 
the NNN spins. In this case of finite systems and with 
chosen values of the exchanges, due to the quantum fluc- 
tuations the values of the order parameters dp and dAp 
deviate from the classical values 0(— 0.75) in the region 
a > Qfc- The oscilations oi dp and dAF (quasi-plateau's) 
at finite N in the region a > are the result of level 
crossing between the ground state and excited states of 
the model. 

To obtain additional insight into the nature of differ- 
ent phases, we have also calculated the x,y, and z com- 
ponents of the dimer-order parameters. We have found 
that dp{d^p) is very close to the saturation value 0.25 
in the region a < The dimerization perpendicular 
to the X axis remains small and close to zero in com- 
plete agreement with the magnetization results. As soon 
as the antiferromagnetic exchange J2 increases from the 
critical value dp{dj^p) jumps to zero. However, all 
components of the F-dimer order parameter (d^^'^) re- 
main close to zero in the region a > Oc- Which shows 
that there is no long-range ferromagnetic order in the 
region a > ac- In contrast, components of the AF- 
dimer order parameter (d^p^) smothly change from al- 
most zero to the value —0.15. Due to the quantum fluc- 
tuations induced by the ferromagnetic exchange Ji, the 
value of these components deviate from the saturation 
value —0.25. 

Thus, our numerical results show that the ground state 
phase diagram of the frustrated ferromagnetic spin-^ 
chain for small values of the antiferromagnetic exchange 
(a < 0.5) contains, besides the gapped ferromagnetic 
phase, the AF-dimer phase. Each phase is characterized 
by its own type of long-range order: the ferromagnetic 
order along x axis in the ferromagnetic phase; the AF- 
dimer order between NNN spins in the AF-dimer phase. 



III. THE SCALING BEHAVIOR OF THE 
GROUND STATE ENERGY 



" 77 XI {~^n-^n+-2j • (7) 
n 

It is clear that the parameter dp{dAF) is the F(AF)- 
dimer order parameter. In Fig. [2K(b) we have plotted 
F(AF)-dimer order parameter as a function of the pa- 
rameter a for the chain with | Ji |= 1.0 and different 
values of the chain lengths N — 20, 26. As is clearly seen 
from this figure, ioi a < a^ dp and dAF are very close 
to 0.25, which confirm that the ground state of the sys- 
tem is in the fully polarized ferromagnetic phase. For 



The finite size scaling method is a way of extracting 
values for critical exponents by observing how measured 
quantities vary as the size L = Na (a is the lattice spac- 
ing and we will consider to be one) of the system studied 
changes. In fact, this method consists of comparing a 
sequence of finite lattices. The finite lattice systems are 
solved exactly, and various quantities can be calculated 
as a function of the lattice size L, for small L. Finally, 
these functions are scaled up to L — > oc^^. Two steps 
are needed before these ideas can be realized. First, one 
needs a procedure for solving the finite lattice systems 
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exactly. Second, one needs a procedure for extrapolating 
from finite to infinite L. In the step one, we have used 
the Lanczos method to obtain the ground state energy. 
We also checked our numerical results by the modified 
Lanczos method"^^. Using the modified Lanczos method 
one can get the excited state energies at the same ac- 
curacy as the ground state one. We did not find any 
irregular size dependence of the ground state energy in 
our numerical results. In the following, we present our 
finite-size scaling approach for the ground state energy. 

Using Lanczos method, we can compute the ground 
state energy as a function of the chain length N and 
the parameter 7 as Eo{N,'y). We have implemented the 
modified Lanczos algorithm on finite size chains (N — 
10, 12, 14, 28) by using periodic boundary conditions 
to calculate the ground state energy as a function of the 
parameter 7. 

In the case of 7 = 0, the spectrum of the ID F-AF 
Ji — J2 model is gapless. The ground state energy in the 
thermodynamic limit behaves as Eq. ^ . By checking the 
behavior of the function Eo{N, 7 = 0) as a function of 
{Eq = AN"), found the best fit to our data yielded A = 
—0.1875 and a — 1.0, which shows very good agreement 
with the analytical result Eq.®. 

Now let us introduce our finite size scaling procedure 
to find the correct critical exponent of the ground state 
energy in the vicinity of the critical point ac = j. First, 
we write the scaling function f{x) as the following ex- 
pression, 

N{Eo{N,^)-Eo{N,^ = 0)) = f{x), (8) 

where x = Nj^ is a scaling parameter. As expected, the 
behavior of this equation in the combined limit 

N — > 00, 7 — > (x » 1) (9) 

is consistent with Eq.(j4]). Thus it can be assumed that 
the asymptotic form of the scaling function is 

fix) ~ xt (10) 

and the ^-exponent in the large- a; regime (a; S> 1) must 
be equal to one — 1). Then we get in the large-x 
regime 

lim f{x)^N{EoiN,-f)~EoiN,0))^x. (11) 

►oo(2;>l) 

This equation shows that the large- a; behavior of the scal- 
ing function ,f{x) is linear in a; = Nj^ where the scaling 
exponent of the ground state energy is /3. We should note 
that in using the Lanczos method we are limited to con- 
sider the maximum value of A'^ = 30*^^. Moreover, since 
the scaling behavior is restricted to the limit 7 — > 0, we 
should consider as soos as possible very small values of 
7 < 0.002. Therefore, the value of a: cannot be increased 
in this method. However, we are not allowed to read the 
scaling exponent of the ground state energy which exists 
in the thermodynamic limit [N — > 00 or a; ^ 1). Thus, 
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FIG. 3: The value of the ground state energy Eo{N, 7) versus 
the parameter 7 close to critical point 7c = 0. The results 
reported are for chain lengths = 20 and best fit is obtained 
by using equation (So(iV,7) oc 7*) with S = 0.999 ± 0.001. 

we have to find the scaling behavior from the small- a; 
regime. According to our numerical computations where 
N < SO, the small- a; regime is equivalent to very small 
values of the parameter 7. In this case the ground state 
energy of the finite size system basically represents the 
perturbative behavior— 

Eo{N,j) = B(°)(A) + B(i)(A)7 

+ B'-^\N)-/^ + ... (12) 

The effect of higher-order terms can be neglected for 
7 < 0.002 to a very good approximation. The first coeffi- 
cient in the perturbation expansion B^-^\N) is the same 
as Eo{N, 7 = 0). To find a relation between other coef- 
ficients and correct critical exponent of the ground state 
energy, it is more convenient rewrite Eq.® as'^® 

EoiN,j) - EoiN,0) ^ g{iN^^f')j), (13) 

where f{x) — Ng{x). This imphes 

— — ° I = Nf X constant, (14) 

where m is the order of the leading term in the pertur- 
bation expansion. Using Eq. p^ we obtain 

B^"'\N)ocNT. (15) 

Now, if we consider the large-A behavior of _b(™) (A^) as 

lim B^'^^N) -aiN", (16) 
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IV. SUMMARY 




E_i I I I I I I I I I I I I I I I I I I I I 

10 15 20 25 30 

N 

FIG. 4: The value of the scahng function B'^^ (iV) versus the 
chain length iV = 10, 12, 14, 28. The best fit is obtained by 
using Eq.([T6| with 6 = 0.62 ± 0.01. 



we find that the critical exponent of the ground state 
energy is related to the ^-exponent as, 



The above arguments suggest that we should look for 
the large- iV behavior of the coefficient B'^"^\N). To do 
this, in the first step we plotted in Fig. [3] the ground 
state energy EQ{N,-i) versus 7 [0.0001 < 7 < 0.0002] for 
a fixed size iV = 20. The best fit to our data is obtained 
with 7 = 0.999 ± 0.001 (£;o(iV,7) oc 7''), which shows 
that the first nonzero correction in the perturbation ex- 
pansion is the first-order (m = 1). We have also im- 
plemented our procedure for different values of the sizes 
N = 10, 12, 14, 28 and found the same results for m as 
we expected. 

In the second step, we fitted the results of the ground 
state energy £'0 [N, 7) to the polynomials for 7 close to 
7 = as Eq. fT^ up to m = 1. Using this procedure we 
found the coefficient of the first-order correction pertur- 
bation, B^^\N), as a function of N . Then we plotted 
in Fig. |4] the function B'^^\N) versus N. The results 
have been plotted for different sizes N — 10, 12, 14, 28 
to derive the ^-exponent defined in Eq. pB]) . We found 
the best fit data for 6 = 0.62 ± 0.01. Therefore, using 
Eq. p?)) we have computed the ground state energy ex- 
ponent (3 = 1.61 ±0.01. Our numerical results show very 
good agreement with the exponent derived in the theo- 
retical point of view, Eq.(|3]). 



To summarize, we have studied the ground state phase 
diagram of the frustrated ferromagnetic spin-1/2 chain 
for small values of parameter a < 0.5. We have imple- 
mented the Lanczos method to obtain the ground state 
energy in small chains. The modified Lanczos method^ 
is used also for checking the numerical results. Using 
the exact diagonalization results, we have calculated the 
various order parameters and spin structure factors as a 
function of the parameter a. It is found that the sponta- 
neous magnetization jumps to zero at the critical value 
a = ac = 1/4. Increasing the antiferromagnetic ex- 
change J2 from critical value ac, the system smoothly 
undergoes a transition from a ferromagnetic phase into a 
phase with the dimer ordering between the NNN spins. 

On the other hand, it is believed that the ground state 
energy behaves as Eq ~ 7'', where /3 is a critical expo- 
nent. From the classical approximations and the spin- 
wave theory, it had been obtained p = 2. Using the vari- 
ational approaches and perturbation theory, it had been 
shown that the quantum fiuctuations definitely change 
the critical exponent and /? = 5/3. On the other hand, 
it had been believed that the exact diagonalization of fi- 
nite chains shows a complicated iregular size dependence 
of the ground state energy, which makes the numerical 
estimation of the critical exponent (3 impossible^^. 

In this paper, we have used the finite-size scaling ap- 
proach to investigate the critical exponent of the ground 
state energy. To estimate the critical exponent of the 
ground state energy, we have introduced a proper scaling 
function f{x) as Eq.®. The scaling variable is defined 
as a; = Nj'^. According to our approach the right scal- 
ing exponent of the ground state energy gives a linear 
behavior of the scaling function f{x) versus x for large 
X. But, the Lanczos numerical results are not able to get 
the large- a; behavior. 

To find the correct critical exponent of the ground state 
energy in the small- cc regime {x <C 1), we have plotted 
the best fit to the data of the scaling function B^™'\N), 
which is the coefficient of the m-order perturbation ex- 
pansion of the ground state energy. The critical exponent 
of the ground state energy is computed with the relation 
between the divergence of the leading term {B^"^\N)) in 
the perturbation expansion and the scaling behavior of 
the ground state energy (Eq. (fT7|) ) . Our numerical results, 
confirm that the quantum fiuctuations are very impor- 
tant and change the critical exponent from the classical 
value and f3 — 1.61 ± 0.02, in good agreement with the 
analytical results (Eq.Q). 
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